Coherent control of lattice deformations in quantum wires by optical self-trapping 
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We investigate theoretically a semiconductor quantum wire under the effect of an intense off- 
resonant cw laser. We show that in the regime of strong exciton-phonon coupling the light-dressed 
ground state of the wire reveals symmetry-breaking features, leading to local lattice deformations. 
Due to the off-resonant nature of the excitation, the deformations are reversible and controllable by 
the intensity and frequency of the laser. We calculate the light-induced strain forces on the lattice 
in the case of an organic quantum wire. 



When an atomic system or a semiconductor is irradi- 
ated by a pump laser in the transparency spectral re- 
gion, it responds to the field with its dynamic polariz- 
ability and gains a polarization energyi In the case of 
semiconductors, this dynamic Stark effect is well under- 
stood in terms of the creation of virtual electron-hole 
pairs (excitons) across the band gap by pump photons, 2 
and it explains the excitonic blue shift experimentally 
observed^ It has been recently pointed out that the dy- 
namic polarization can be used as a quantum control tool 
in semiconductors. For instace, the virtual excitons cre- 
ated by an off-resonant pump field interact with spins 
localized by impurities or quantum dots and can create 
a local magnetic field^ induce spin-spin coupling, 5 ** and 
paramagnetic to ferromagnetic transitions^ 

In this paper we will show that the dynamic optical po- 
larization can be used to induce strain forces and control 
lattice deformations in a semiconductor quantum wire. 
If the virtual excitons created by the pump field interact 
strongly with phonons in the lattice, a self-trapping* of 
the optical polarization results. This induces local forces 
on the lattice. As in the case of the optical spin control, 
there is no absorption of energy since the pump photons 
are tuned in the transparency region: the effect is due to 
radiative (stimulated) corrections to the ground state of 
the system. The coherent nature of the effect makes it 
reversible and finely controllable with lasers. In principle 
many different materials have the exciton-phonon cou- 
pling strong enough for this optical self-trapping effect. 9 
Here we focus on an organic one dimensional system, 
polydiacctylcne (PDA) j** which has been studied for its 
strong nonresonant optical nonlinearities and exciton- 
phonon effects. Phonon-mediated optical nonlinearities 
have been observed in this material. 11 

Recently-developed experimental techniques can de- 
tect light-induced lattice displacement in molecules and 
semiconductor systems. Examples include pump-probe 
electron diffraction** and ultrafast x-ray absorption spec- 
troscopy!** These techniques can be extremely sensitive 
and measure laser-induced lattice dynamics within pi- 
cosecond and milli- Angstrom resolution. 14 PDA has been 
already identified as a good system where the vibrational 
dynamics after optical excitation could be observed^ 
and the possibility of addressing optically a single poly- 
mer chain of this material has been recently demon- 
strated]** We are focusing here on the steady state regime 



of the polymer driven by a cw or nanosecond laser. This 
should be easier to address experimentally, yet it contains 
rich and unexplored features related to light-matter in- 
teraction in strong coupling. We remark that, in contrast 
to more common schemes based on photothermal effects, 
the coherence of the electronic excitation is a key element 
in this light-induced effect i** 

The coherent many-body ground state of the wire in 
the presence of the light field can be parametrized as a 
BCS-like wavefunction that depends on variational pa- 
rameters. By a functional minimization of the ground 
state energy with respect to these variational parame- 
ters, we obtain a non-homogeneous nonlinear equation 
that describes the distribution of the optical polarization 
in the wire. Unlike the case of a single exciton-polaron, 18 
the total polarization is not fixed but is determined by the 
intensity and frequency of the laser field. The nonlinear 
equation is solved numerically to calculate the distribu- 
tion of the polarization and therefore the strain forces 
acting on the wire. Recent esperiments on single PDA 
chains 16, 19-20 suggest that free excitonic states have a 
Bohr radius of 1-2 nm, corresponding to 2-4 monomer 
molecules, and a mass of 0.3mo. We use these param- 
eters to describe the wire with a tight-binding model. 
The model takes explicitly into account the exciton-light 
coupling. The hamiltonian can be written as 

H = H SSH + H X l ; (1) 

where the first term has a Su-Schrieffcr-Hccgcr (SSH) 21 
form 



H 



SSH 



E^+E^( u «+ i - M ») 



c , 



- ^ tn+l,n{B n+l B n + E^Bn+l) ) (2) 
n 

and the exciton-light interaction is described by 
Hxl = ^(Bl+B n )+Y,(E g +2t -hu; L )BlB n . (3) 

n n 

M , C, u n and p n are the mass, elastic constant, to- 
tal displacement and momentum of the n-th site of the 
tight-binding chain. _B+, B n are operators of creation, 
annihilation of excitons in a singlet spin state. The hop- 
ping term is i„+i,„ = t - j(u n+ i - u n ), where 7 is re- 
lated exciton-phonon deformation potential interaction: 
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D = 2ja. a is the site separation in the tight-binding 
chain taken equal to the excitonic Bohr radius (1.5 nm for 
PDA) . The value of to is determined by the exciton effec- 
tive mass m as to — h 2 /2ma 2 . We treat the exciton-light 
coupling (Hxl) scmiclassically, and the rotating wave 
approximation is assumed. The shift 2to is introduced 
in such a way to define E g as the k = exciton energy 
at the bottom of the band, and lul is the frequency of 
the laser, = d£o is the Rabi energy due to the dipolar 
coupling of excitons with a light field of amplitude So- 
We remark that PDA has a large exciton binding energy 
(of the order of 500 meV) and also shows vibronic reso- 
nances spectrally localized at about 200 meV below the 
free exciton peak3 If the optical detuning S = E g — Tilul 
is small compared to these quantities, both free e-h pairs 
and vibronic excitations can be neglected in the model. 

In contrast to the polyacetylene case^l the ground 
state of polydiacetylene is a non-degenerate state that 
we can express in the form 



l$n 
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(4) 



meaning that the wire has no excitons. However, in the 
presence of light, the coherent light-dressed ground state 
will change, and we write it using the ansata^ 



|*n) =Y[(a n \0}n + l3ne^\l)n) 



(5) 



where a ni /3„, and (p n are variational parameters, and 
a n and f3 n are subject to the single-occupancy constraint 
\ot n \ 2 + \Pn\ 2 — 1- |0)„ and |l) n denote the number of 
excitons in the site n. Due to the hopping, the variational 
coefficients depend on the index n. We introduce the 
phases tp n to make a n and (3 n real. The mimimum of 
the total energy occurs when the phases are constant, 
making all the two-level oscillators in the chain mutually 
coherent^ 3 . Using the wavefunction in Eq.© and the full 
Hamiltonian in Eq. ©, we calculate the total energy. By 
introducing the optical polarization function 



ip n = -(B n ) = 2a n /3 n , 
the energy can be written as 
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(terms that do not depend on tp n have been dropped) 
where we have defined 6' = 6 + 2t - We also have set 
tp n = 7T, and assumed a n /3 n > and \a n \ > \/3 n \ ■ By 
functional derivation with respect to ip n we obtain 



-t 



t n ,n-\i>n-\ —fl + S' 



0, (8) 



which in the continuum limit reads 



-to [Wn + 27^^ - n + s -j=== = o. 

(9) 

In order to have a closed equation for ip n wc need to find 
du n /dn, which depends self-consistently on the distri- 
bution of the polarization along the chain in the steady 
state. This can be calculated using the same approch of 
Ref. by writing the classical equations of motions for 
the wire 



Mu n = 



d(H) 



C(u n+ l - U n ) - C{u n - Un-l) + F n 



( 10 ) 

In the continuum limit the strain force on the site n, F ni 
is related to the gradient of the polarization density as 



F n = 



7 9\lp r , 



2 dn 



(11) 



Since we are considering a stationary solution in the 
presence of the light we set Eq. H10JI to zero, which 
gives the local displacement due to the strain force as 
u„ = -C- 1 J J F n d 2 n. Then from Eq. (TTJ we find 

£~i|*.|« + < * (12) 

where A is a dimensionless constant of integration that 
gives finite-size effects. We choose A = 0, which implies 
that the total lenght of the wire is not fixed. This choice 
does not affect the results in the limit of a wire much 
longer than the extension of the self-trapping region. The 
force constant C can be expressed in terms of the sound 
velocity S as C = S 2 M/a 2 where S = 2.5 x 10 3 m s" 1 
for PDAiS Finally, from Eqs. © and we get the 
equation for ip n 



/ 1 \ d 2 ip n D 2 o , ip n 



I 

Notice that we are dealing with a system where the quan- tion density in the wire, i.e. the number of virtual ex- 
tity tp n is not normalized. The amount of total polariza- citons, is determined by the external field parameters f2 



3 



and 5. The first term in Eq. (|13f) describes the blue shift 
of the exciton (Pauli blocking). The cubic term, related 
to the exciton-phonon deformation potential interaction 
D, is responsible for the self-trapping of the polarization. 
It corresponds to an attractive potential proportional to 
the polarization density. The last term takes into account 
the effect of the finite optical detuning of the laser field. 
In the absence of light (i.e. for £1 = 0), the solution to 
Eq. ifTsTl is ip n — corresponding to the trivial solution 
in Eq. |0J. In the absence of hopping between the sites 
(<o = 0, D = 0), we have spatially separated excitations 
with no constraints applied^ The solution in this case 
corresponds to the one of the dynamic Stark effect in 
atomai 



At resonance (<5 = 0) Eq. i|14fl gives the maximum occu- 
pation ip n = 1, corresponding to half exciton per site, and 
the light-dressed ground state is 2~ 1 / 2 (|0) — |1)) on each 
site. In the limit of small excitation ifi n <§; 1, Eq. (|13|) is 
similar to a one-dimensional nonlinear Schrodinger equa- 
tion with a cubic attractive nonlincarity. It is known that 
a Bose-Einstein Condensate (BEC) with an attractive 
nonlinearity breakes symmetry in ID so that the ground 
state does not correspond to a constant amplitude so- 
lution^ In our case, these real symmetry-breaking so- 
lutions for the ground state describe the optical self- 
trapping effect. In the small excitation limit some an- 
alytical forms for the solution of Eq. 1)1 -iH can be written 
in terms of Elliptic functions using Cayley's transforma- 
tions^ It is interesting to remark that the Hamiltonian 
in Eq. can be mapped to an XY model for spins cou- 
pled to phonons. The Rabi energy describes an external 
magnetic field in the XY plane while the optical detuning 
corresponds to a magnetic field in the Z direction. An 
effective nonlinear interaction between spins analogous 
to our cubic term was found in the case of a classical 
Heisenberg spin chain coupled to phonons^ 

In order to study on the same ground the inter-cell 
hopping and the presence of electromagnetic field we 
solve numerically Eq. i|13|l • This allow us to calculate the 
forces on the lattice for an arbitrary value of the external 
field, provided the system remains coherent. To solve the 
equation numerically, we use the steepest descent method 
of functional minimization^ This method is efficient in 
solving numerically Gross-Pitaevskii equations for BECs, 
which are in the same class of Eq. Ijl3|l , 27 We start from 
an initial trial function ?/Vi(0), then ipn(r) is evaluated in 
terms of the equation d ^£^ = —Titp n (T) , by choosing an 
arbitrary step At in the fictitious time variable r and by 
iterating according to ^ n {T + At) VVi(t) — ArTiipni''') ■ 
The nonlinear operator TL is defined in such a way that 
Hipnij) — corresponds to Eq. ljT3|l. At convergence, 
the condition ^f - = gives the solution to Eq. (|T3*|l . 

The results of the numerical solution for tp n are shown 
in Fig. ^ The linear density of the wire is 1.4 x 10~ 14 
g/cmii&. The sum of the deformation potentials for the 
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FIG. 1: (Color online) Polarization density \ip n \ 2 for 8 /to = 
0.05, Q/to = 10" 4 (green line), Q/t = 1.99 x 10 -3 (red line), 
and Q/t = 2.06 x 10" 3 (blue line) 
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FIG. 2: (Color online) Force as a function of the site number 
for 8/to = 0.05, Q/to = 10" 4 (green line), Q/t = 1.99 x 10~ 3 
(red line), and Q/to = 2.06 x 10" 3 (blue line) 



conduction and valence band is D = D c + D v = 6.1 eVi^ 
The dipole momentum of PDA is 10.1 eAj2£ and the en- 
ergy of the excitonic (singlet) state is 2.282 eVii&The nu- 
merical solution is obtained in a chain with 100 sites and 
periodic boundary conditions are used. The total lenght 
of the chain is smaller than the wavelength of the light 
and we can neglect the effect of a finite photon wavevec- 
tor. The experimental value of the excitonic linewidth at 
10 K is 350 /ieV,— which is much smaller than the values 
for the Rabi energy and detuning we use in the calcula- 
tions. We have also checked that the exciton-polaron 
binding energy is always smaller than 5. Fig. ^ shows 
\ip n \ 2 for three different values of the Rabi energy and 
the same optical detuning S/to = 0.05 (to = 70 meV) as 
a function of the site number. We estimate that the value 
of the Rabi energy needed for self-trapping corresponds 
to laser intensities of about 2.5 kW/cm 2 , which is rea- 
sonable for nanosecond pulses. At the lowest value of the 
Rabi energy the polarization is small and homogeneously 
distributed along the wire. For fl ~ 150fieV, the solu- 
tion is localized due to the phonon-assisted self-trapping. 
When the Rabi energy increases slightly, the polariza- 
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tion distribution becomes broader and flatter at the top. 
This saturation is related to the intrinsic fermionic na- 
ture of the excitons. Notice that the total polarization 
increases when the Rabi energy increases. A gradient in 
the density of polarization produces a force on the ions 
according to Eq. The force is stronger at edges of 

the |VVi| 2 distribution as seen in Fig. [21 and is positive 
to the right from the center of the chain and negative to 
the left. The quantity that is measured in electron or 
x-ray diffraction experiments would be du n /dn, i.e. the 
local deformation due to the light-induced changes in the 
lattice. For the parameters considered above, this quan- 
tity is in the range 10-100 mA, which is reasonable for 
an experimental observation of the effect. 

We have neglected exciton-exciton Coulomb effects 
which may affect the picture given here, especially in 
the high excitation regime. This deserves further in- 
vestigations. However, it is known that the BCS vari- 
ational ansatz in Eq. JSJ provides a good description 
of the ground state of an electron-hole system both in 
the low excitation regime (excitonic BEC) and in the 
high excitation regime (BCS pairing) m& The fact that 
the electron-hole system here is driven by an external 
coherent field (laser) gives an additional justification to 



this ansatz since light-dressing can reduce Coulomb fea- 
tures beyond mean field like screening and broadening^ 
Finally, exciton-exciton interaction in the dynamic Stark 
shift can be neglected when the optical detuning is big- 
ger than the biexciton binding energy^ therefore its effect 
can always be reduced by increasing the optical detuning. 

In conclusion, we have investigated the effect of an in- 
tense off-resonant laser field in polydiacetylene chains. 
We have extended the SSH model to include the effect of 
a control laser field semiclassically. We have obtained a 
nonlinear inhomogencous equation describing the spatial 
distribution of the optical polarization in the chain in the 
steady state. We find solutions describing self-trapping of 
the polarization which saturate when the intensity of the 
field is increased. The results suggest that local steady 
state lattice deformations can be finely controlled by the 
intensity and frequency of an intense laser field. This con- 
trol scheme could have interesting applications to other 
organic materials where a strong exciton-phonon interac- 
tion exists, like e.g. in J-aggregates, or DNAi^I 
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